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1. INTRODUCTION 

The performance prediction for multicarrier operation of a nonlinear 
satellite repeater is essential for the optimization of a telecommunications 
system which operates in the multiple access mode with frequency distribu- 
tion. 


The spectra of the signals running through the equipment undergo 
deformations due to intermodulation in the. progressive wave amplifiers 
(TOP), low noise amplifiers, mixers or filters. Depending on their mode 
or operation point, this equipment behaves more or less like nonlinear ele- 
ments. This produces amplitude distortions and phase distortions which are 
detrimental to the signal being transmitted. Basically, this equipment has 
two kinds of nonlinearities [1]: 

— a phase nonlinearity; 

— an amplitude nonlinearity. 

In most cases these nonlinearities have no memory p can be modeled, and 
can be measured: 


Numbers in the margin indicate pagination of original foreign text. 


- They have no memory because the output voltage of the element 

only depends on its direct gain and the input voltage: this is 

the case of the TOP and the filters; 

- They can be modeled because they can be analytically simulated 
by complex polynomials; 


— They can be measured using a network analyzer in the case where 
the number of carriers is not too large. 


Schematically, in the case of an input signal consisting of two 
carriers of equal amplitude and separated by AF, the problem can be formu- 
lated in the -Following way: 


{? ° 


- 2 carriers 


4 * 



t . iA T 


-L 


_ ___ . 

HA 


1 


_ 2 carriers and 
intermodulation 
products on the 
. order of 3, 5, . . . 


/4 


In the present report we will describe a "source program" which allows 
one to calculate the following from subprograms of the COMLIB [2]: 

— the complex envelope of the output signal; 

— its spectrum; 

— the carrier ratios with intermodulation products; 

— the carrier ratios with intermodulation noise in a frequency band 
(C/N). 
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The input signal can consist of as many lines as necessary. 


The results of simulation will be compared to those obtained in ex- 
periments. A direct application consists of a low-noise amplifier which 
operates in the frequency range between 3.7 to 4.2 GHz [3], which is de- 
signed for a 6-4 GHz repeater. The input signals consist of two or three 
carriers and of a noise density associated with a single carrier of low 
amplitude. The simulation program will also be applied to the case of an 
input amplifier for a direct TV satellite, for which the transmission/re- 
ception coupling introduces parasitic signals at the input which inter- 
modulate with the useful signals. Such a program, therefore, is a working 
tool which is useful for "conceiving systems" and is also useful to the 
"circuit engineer" because, from overall system spec? fications (for ex- 
ample, minimum C/1) , it can help in the selection of the component or the 
equipment and can also specify its optimum operational point. 


2. STARTING HYPOTHESES REGARDING THE 3.7 - 4.2 GHz AMPLIFIER 

The nonlinear model of the amplifier considered consists only of 
amplitude and phase nonlinearities, In particular, we will ignore the 
amplitude/frequency nonlinearities which would lead to a numerical resolu- 
tion system which would be too complicated. 


The amplifier is considered as having no memory. That is to say, 
there are no reactions at the output to the input. In practice, there is 
a complex inverse gain S^ but in the first approximation we can assume 
that it is sufficiently small compared with the direct gain S2^, so that 
it can be ignored. In addition, the amplifier is assumed to have a limited 
band so as to neglect the influence of the input signal harmonics. 


The following is the equivalent diagram for the amplifier: 


passband fil- amplitude and passband fil- 
ter phase non- ter 

linearities 


3. THEORETICAL ANALYSIS 

For more clarity, we seem to prefer to collect several terms about 
complex signals and complex gains developed by M. Sombrin in Report No. 
77/026/CT/TCS/LT [1] . 


3.1 Concept of Complex Gain 

The amplitude and phase nonlinearities have the following appearance: 




1 dB point for gain compression 
where ‘ : input signal amplitude 

: out P ut signal amplitude 
^ : phase difference between output and input 
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We can then introduce a complex gain of the amplifier which takes into 


account these nonlinearities: 
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In the case of a low noise amplifier, we arbitrarily normalize the 
complex gain to 1 for the 1 dB compression point: 

This normalization allows one to generalize the method and to apply 
any kind of nonlinear element: mixer, limiter, assuming that the coordi- 
nates of the 1 dB compression point are known. 


3.2. Representation of Signals as a Function of Time 


n 


The numerical representation of a signal as a function of time re- 
quires the sampling of this signal. The number of samplings will become 
prohibitive at several gigahertz and because of the theorem of Shannon. 


It can be shown [1] that the signal can be represented only in the 
form of its complex envelope, that is if S(t) is a real signal which is 
modulated in amplitude and phase: 


S(t) s ECtJ.coo [u>*fc +• tb (fc)l 


set) 2 %i j 


The complex envelope Ec(t) around fo is: 

k'O. © ,>c<0 = X(k) + I VO) 


We then have: 


s(t). = i yetj] . e 1,A,t J 
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If B Is the passband of the signal, the maximum frequency of Ec(t) 
is then B/2. The sampling will be done at the frequency B Instead of at 
fo. The calculation time is then reduced accordingly. 


3.3. More General Case 


If we apply the nonlinearity over the complex envelope of the signal, 
we then find: 


— real signal: eft) » E{t) eosQu> 0 fc +• 

— complex signal around Se(t) * . a 

— complex gain of the nonlinearity: ° 


JVC?.) 


— complex envelope at output: 


. Sc(u; » E(t) . c 


E(t).S [£^j • C 


(«?•>) 3 


real output signal: 
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from measurements. 


3.4 . Properties of Fast Fourier Transform 

The transfer from the complex envelope of the signal to its spectrum 
and inverse operation implies the use of a Fourier transform. In order to 
reduce the calculation time and increase the accuracy of the results, at the 
present time numerical calculus methods use the Fast Fourier Transform method 
(FFT) . 


3.4,1. Time-Frequency Correspondence for FFT 


The two diagrams given below give an example of using discrete FFT for 
calculating discrete spectra and the complex envelope of the signal; 
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These time-frequency relationships, therefore, satisfy the following 
relationships : 
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4. NUMERICAL SOLUTION PROCESS 
4.1. Synopsis 

Below we will give an organization diagram of the subprograms which 
are called by the main program LINA and their mathematical significance. 
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4 . 2 Description of the Main Program LINA 

LINA carries out the operations described above If the user has supplied 
a certain number of data about the following: 

— the number of samples: N; 

— the input spectrum (phase amplitude and frequency of each line) ; 

— the polynomials which model the nonlinearities; 

— the coordinates of the 1 dB compression point; 

— the input power range; 

• — at the output, the various plotting options desired (either as list- 
ings or as Benson prints) . 

Appendix 1 gives the details about the data cards and their input format, 
as well as the FORTRAN of the subprograms which do not belong either to COMLIB 
or to CNESLIB. 


4 • 3 Precautions ff * s. Use 

A certain number of precautions must be maintained In order to obtain 
correct results. First of all, the number of samplings must be sufficient 
so as to contain the lines of the inlet spectrum and the lines which are 
suspected In the output spectrum. 

An excessively small number of samplings can lead to a folding 
phenomena of the spectrum which is related to the periodicity of the dis- 
crete Fourier transform. This phenomenon is illustrated by the following 
figures: 


Normal spectrum, 
N Is sufficient 


Incorrect output spectrum, 
M is insufficient 



In [1], the folded-back zone of the spectrum is added to the real 
spectrum. At the end of the spectrum there is then a recombination in ampli- 
tude and phase of these lines, which can lead to substantial errors. 


In practice, the number of samples must cover twice the pas sb and of the 
input signal. In addition, the number of samples must be a power of 2: 


5. 


RESULTS 


5.1. 


Case Where the Input Signal Consists of 2 Carriers of Equal Amplitude 


5.1.1 Experiment s 


Two carriers spaced by 5 MHz around 4 GHz are reentered In the end 
stage of an amplifier with an Intermediate frequency, operating In the fre- 
quency bond between 3.7 to 4.2 GHz. 


The intermodulation products created by the phase and amplitude non 
linearities have the following form: 

2. - t rn£ ± r 

/<• h - U* 


where: £1 and f 2 : the two carriers 

m * 0, 1, 2... order of Intermodulation products 
n - 0, 1, 2... " " " 

If we ignore products with an order greater than 7, we obtain an output 
spectrum having the following form: 



The test equipment used Is the one described In Figure 1. 
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Figure 1 


The ratios C/I^ are measured for different values of total input power 
and are plotted in the curves of Figure 2. 

The test complex allows one to determine the amplitude-phase relation- 
ship of the amplifier, and was previously described in Report No. 76/061/DSES/ 
RT/TTI/M0/E2 . The amplitude and phase curves for this amplifier are shown in 
Figure 3. 


5.1.2 Simulation 


We select 256 samples for representing the input signal. The output 
spectrum is given in Figure 4 for a total input power of - 20 dBm. 

The intermodulation products calculated according to the previously de- 
scribed method were plotted in Figure 2 and compared with the measured values. 
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Figure 2 
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Figure 3 
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?igure 4. Output spectrum of the amplifier 
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A first analysis of these curves shows a deviation between the measured /17 

and calculated values . As far as the "low level" points (P^ n > -10 dBm) are 
concerned, the mathematical model selected is found to be insufficient. In 
effect, the amplifier in this range behaves like a memory element, a hypothesis 
which has been Ignored from the start. In effect, this is not objectionable 
because the input power range in which we are interested is lower. 

Also, points measured for C/I greater than 60 dB have errors because of 
the fact that \?e do not know the dynamics of the spectrum analyzer used. 


Nevertheless , we find the general variation of the theoretical curve 
with a slope of 2, which results from the ratio of a curve having the slope 
1 (Ps = f (Pin) and a curve having a slope 3 (Px^ = f(Pin)) . 

In figure 5 x?e show the curves calculated from the program for different 
intermodulation products of order 3, 5, 7 and 9 as a function of the total in- 
put power. 


We can note that the interception point is the same for all of them and 
that its coordinates are: 


f *in M 5 dFm 

Interception point J 

l Pr» « 18 dBm 


Products of order greater than 3 are negligible, which justifies that 
we only consider products on the order of 3 in this study. 

for the hypothesis of a ratio of C/I of 50 dB with 2 carriers equivalent 
to 28 dB in the favored channel among the 12 equidistant channels with identi- 
cal levels, we obtain the maximum admissible power of the component utilized 
(here, a bipolar transistor CG126) or, in other words: 

Ps - t- 10 dBm 
max 
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5.2 Input Signal Consisting of 3 Carriers of Equal Amplitude 
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The Input spectrum has the following appearance: 
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This configuration of lines allows one to dissociate products on the order 
of 3 and to add their power levels in the calculation of the global C/l for 
these products. The following figure shows an output spectrum obtained under 
these conditions if we are limited to products on the order of 3. 

hr 


■m- 


t 


* 

e*r- 


I 


I 








1 

‘Zn 

o-ji 

l 

1 


< vT 
1 

ala 



f 


5.2.1 Experiments 


The test equipment used for this experiment is similar to the one de- 
scribed before. In addition, the wave generators must have exceptional 
stability in terms of frequency. The frequency separation (f2 - fl) = AF 
was selected equal to 5 MHz in order to provide good discrimination of the 
rays for the analyzer, without excessively widening the frequency spectrum. 
Nevertheless, the total frequency band is still sufficiently small so that 
one can neglect the gain-frequency nonlinearity. 
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5.2.2 Simulation 


The three lines have the same amplitude but a phase which is equally 
distributed over the Interval (0.2 tt). Figure 6 shows the spectrum cal- 
culated for an input power of -25 dBm. 

Figure 7 can be used to compare the calculation results and the ex- 
perimental results for total input powers of -22 dBm and -13 dBm. 

We can see that these results are quite close aiad, as predicted, the 
largest product is the one of order 3 having the form (+ fl + f2 + f3) . 

Figure 8 shows the curve (C/I) * f(Pin), where X is the total power 
of the intermodulation products located between f2 and f3. This allows one 
to evaluate the ratio G/X for the hypothesis where the difference between 
f2 and f3 is only 2 AF; this is the least favorable case. 


5,3 Case Where the Input Signal is a Noise Density 

The multicarrier behavior of a nonlinear element can be derived from 
the behavior obtained in the presence of white noise which is limited within 
the frequency band of interest. 

The nonlinearity is then expressed as a ratio of noise powers. Tills 
amounts to finding the ratio of the noise density in a certain frequency 
band Af and the noise density in a window with the same width. The follow*' 
ing figure describes this operation: 
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figure 6 
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Figure 8 
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The carrier ratio for Intcrmodulatlon products In this case Is then*. 


m. 



?» 


¥ 




This method requires the use of a very selective filter of the stop-band typo 
In order to obtain a narrow window with steep edges. Such a filter will re- 
quire a special study. 

We have preferred to associate a sufficiently weak carrier with the 
noise density, compared with the total noise power around which the Inter- 
modulation noise is located. 


Input Output 



5.3.1 Experiments 

Elliptical filters [4] centered around 3960 KHz deliver a noise band 
of 36 MHz when the carrier Is located at 4085 KHz, 

The test equipment is shown in Figure 9. 
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Amplifier-' 



Figure 9 

The measurements were made In a 100 KHz analysis band, a value usually 
used for this type of measurement, This offers an acceptable compromise be- 
tween the rapidity of acquiring a measurement and the minimum function dy- 
namics of the video band of the spectrum analyzer. 

The measurement results are shown In Figure 10. 

5.3.2 Simulation 

In our program, the noise density is simulated by 15 carriers of equal 
amplitude and having random phases. In order to reduce the number of samplings 
and therefore the calculation time, the frequency sampling step was selected 
equal to 1 MHz. 

Figure 11 shows the spectrum obtained at the output for a total Input 
power of -20 dBm. Figure 10 shows- the results of (C/N) measurements and the 
(C/N) calculated values, but refer to a 100 KHz band. 

The deviation between these two curves for the one hand is caused by 
the noise model used (15 carriers with a random phase do not represent a suf- 
ficiently accurate representation of the noise). Also, it is due to the fact 
that the noise power, measurement of the spectrum analyzer has errors due to 
the fact that the exact noise bands of the receiver are not known. 






* 
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Figure 11 
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5 . 4 Application to the Case, of a Low Noise Amplifier Which / 28 

Processes TV Signals 

At the input of a direct TV repeater there are useful signals to be 
treated, as well as parasite signals which come from the coupling of the 
transmitter. Even they are strongly attenuated by Input filtering elements, 
these signals retain an amplitude which is sufficiently large to produce Inter- 
modulation products at the Input amplitude level. 

Simplified, the Input spectrum can be represented In the following way: 


10 useful carriers 


5 parasite carriers from the 
** , s transmitter 




H-KM* 




•9 ftp 


C.SS Mil.- 


4 4 ‘ i? 3 -to 






m a 




Among the 10 first carriers, 5 belong to the French network and the 
others belong to the Luxembourg network. 


We are interested in the spectrum at the output of the input amplifier, 
as well as the (C/N) around each useful carrier which has been shifted before- 
hand . 


The amplifier characteristics were selected similar to those of a para- 
metric amplifier or an FET amplifier, which has the following operating point 
for a C/I of 61 dB for 2 carriers : 

Pin = -59.1 dBm 
PS = -44.35 dBm 
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From these values, the 1 dB compression point calculated is then the 
following: 


rin comp u ' 38 « i dB * 


Po comp “" 24 » 35 dB ® 


A simulation was made using 4096 samples, random phases for 15 carriers 
and a frequency sampling step of 1 MHz. 

Figure 12 shows the ratio (C/I) of the amplifier for two carriers and 
the (C/N) for 15 carriers around carrier No. 7, which is the least favorable. 




29 




6 . 


CONCLUSION 
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The calculation program described within a study of nonllnearltles In 
transistor amplifiers allows one to determine the output spectrum and the 
carrier ratios with Intermodulation products (C/I) in the case where the 
signals are known In the form of their line spectra. 

This program can then be used to evaluate nonlinear elements In the 
amplitude-phase sense, such as TOP, FET amplifiers, filters, mixers, for 
which one ignores the multicarrier behavior or when there is noise. In 
addition, if necessary, one can discriminate between the two kinds of non- 
linearities, amplitude and phase. 

The tests carried out with an amplifier having an Intermediate fre- 
quency of a 6-4GHz repeater showed, compared with simulation results, that 
the selected nonlinearity model agrees especially well with low input power 
levels . 

An additional study carried out jointly .with the Mathematics Division 
will allow a refinement of this model in the case of a nonlinearity with 
memory. 

Finally, application to a case of an input amplifier used in a direct 
TV transponder showed the influence of the parasite signals due to the emis- 
sion on the C/I ratio. 


The author wishes to thank J. Sombrin and H. Donat for their discussions 
and the formulation of the problem as well as its programming. 
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APPENDIX 1 


/3 3 


Meaning of elements contained In the data cards of 
the nonlinearity simulation program LINA. 


Data cards 


Number 

Contents 

Meaning 

1 

N, NP, AIN, ASO, DEL, 
FREQ, FO, Nl, N2, 

PI Comp, P2 Comp 

N: number of samples 
NPi number of carriers 
AIN: total inlet power min, (dBm) 

ASO: total inlet power max. (dBm) 

DEL: Increment in input power. These three 
latter items are referred to the cal- 
culation and plotting of the. output 
spectrum. 

FREQ: total frequency band considered (for the 
N samplings) (MHz) 

^0: central frequency of sampling (MHz) 

Nl) output spectrum, will be plotted as a 
) listing of the sampling 
N2) of the sampling Nl with [illegible[ N2 
PI Comp: input power at 1 dB of compression (dBm) 
PS Comp: output power at 1 dB of comprssion (dBm) 

2 

NIP, AMP, PHI 

i 

NIP: subscript of the first carrier 
AMP: amplitude of the carrier (in volts) 
PHI: phase of the carrier (radians) 

3 

NIP, AMP, PHI 

i. 

Other carriers 

NP + 2 

1 

A(l), A(2) , A(3), 
A(4) 

A (4) . . .A (4).: coefficients of the curve 

output amplifier = f (input amp.) 

NP + 3 

A(5) , . . . A(9) 

A(5) . . .A(9) : continuation of preceding table 

NP + 4 

B(l) , . . . B(4) 

B(1)...B(4): coefficient of the curve 

phase = g (input amplifier) 

NP + 5 

B(5),... B(8) 

B(5)...B(8): continuation of preceding table 
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«r 


Precautions about use : 

Note 1 ; We will give a maximum amplitude of 1 to the carrier 
numbered (2 No * A + 1) 

Note 2 ; Several data packages can be entered once, For each 

data package, a I illegible] card is needed; at the end 
of the last card, a card FIN is required: 


/•End 

I ! ff 




*# *»**;.#* 




2nd data package- 


f (_ blank card 

data package 

/^blank card 


Note 3 : The variables (givens) start with N and must have the 

following format 1^ : 

Example: 12, 30, 25 and not : 12., 20., 30. 

All the other cards must have the free format. 
Example: 2,5, 0.25 E +• 01, etc. 
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/_35 


s? 


PCOG"A h l:M (iN c U,HTnUT,rAPE>;*INrUT, TAPCCbCUPUI ) 

£TUSu CCS N0N-LINEA3I7E5 0, UN AFFIX PS A 

T s A*i$ZSTC c S 8 XFC!.A: v .5S 

cc 4 m cn ’UfSa 

CCMBJ.EX Z(i 3 ;o ,Z 2 Ul-jt) 

C I.M £h*3 Z3f» m3S'J> , Y( 4 155 ) 

cimsvjsos a(?)»3iai» :s\(i.Mi.),£(neA) 

CDi-NCZCN Xf|CP(2?)»<<MM21J 

c m e.ns * a s =u i ( g j n , c sl ?: < in > 

CC m PlE X GZtGC 
INTEGER FIN 
ICA5C*5 
IF*'UNT=S 
CAT A FIM/lCHFIfs 
CCA 3 t 
F CP, F A i 


ca«o* 33 ) :rr-z'i 
(a id ) 


IP( IFF I*|, EC. FIN) 
FI-3.1A1 ‘9 
FACT=1.2 
AY= ] . 


STOP 


S5= 3*0 

sea*:, ■: 

IA=57«£21 
CALL 9 AN SET ( I A ) 

REAC * rN ,NC, tlb , AS3, C5 l, FF- 


i#. k 


» 


- A m t A t*A <* f j c 

*♦ w ^ | u cu K ki; 


, 3,55 


u =2 Y 
NN=N/4 
CC 1 1=1, N 
22 ( Z>b (3 ,i 

cc.’irms 

CC 2 l3l,NP 

READ '‘,NIP,Af'F,P» 4 I 

ALEAsCiANF( JA) 

FH“= (Phj>g*PJ) * A L S A 
VJ P I T 5 * 1 , NIF, A N P ,P(-I 

52 a flvo .r;|K) f ci |, j ) 

tl S3 = AH c * CC S ( FH I ) 

Z2( S :P )sC“PU< (2153 ,33 ) 
CCITIHi: 

HEA3 *i (ill) ,1=1,4) 


3,F*2.M,H2» 

, N i , ') g 


*, t r 


PICCFP (FSCCMp 
,F IC Cii * , FSCONP 




WRITE *, (A (I)*: 


• « » 


4) 


sm D ** , (S( J) , j 


BEAC 

w*:ts 

= c 

W •? * j c 
ft z A C 

^*:ts 

FCAO “ ,1 TF A1 ,I7F4g T5;3 
WRITE *, IT -“A 1, ITAA2, IT *4 
F£AC * » A F I A » il'AX.STr. c ,hC 
►SITE *, AM If) ,AWfl* ( ST HP ,NCI, N 

ifcnci.e:.'!) co rc 121 
nci r •••’(€ 1 


'»(A(I> ,1 = 5,5) 
•, ! A It >,1=5,35 
1 ,«) 

fc , <e u ) ,j = i ,u) 

,(£(«) , JaS.fi) 

-, IE (.J ) , J = 5 , 3) 


,ITPA 4 ,nPA 5 fZ 
2 ,ITFA 4 , J 75 AS, 
l ,M 'i y'(4 i N 5 
i N 4 , K 5 


r A AG 
ITS. A 6 


2 c (NCI 
CC 4 : U B J 


<M- 1 > !' C 

no::: 
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«-> o «“> 


PSA 3 * 

k=:ts ♦ tWPrf-c" 0 
*1 1 :NUM(J)«hUP 

KNU.*l(J)»NU. wo 
J W J * J 

<0 cwriMis 

J 3 1 CCMTms 

CAlCLL ENVELOPE COMPLEX® DU s:G>AL ENTSSE 

VV» ( (PJCC*P-3C, )/??. ) ‘ALCCU3. ) 

Wal./cVPCiV) 

VCs ( ("SCCMP-33. »/2«, l*/LCG( 1C* > 

VC«E>P (VC) 

p»a 

:c?t*3 

cP«ppec/N' 

T*1 » /OF 

IP(ne.Al.cC*l#0"«ITS.Ag.SC.L»CR«tTC43,EC.l) CALL FLCTCCtD 
fi'ut USC-AIM /CEL* 1 1 1 
JJnIFIX 4Y> 

AN s3.3 
CCPps.1. 

FPsPFrtCO (N ,ZZ) 

WsiT6(6»t») «P 
13 fCP.«ATliCV,*FLSS. TOT, IMTlALt(C? M )«J* »F 11 .4 ) 

CCFF*AIN-PP 

cc 

CC 11 s!«l,h 

2 (J ) »Z 2 < J) 

Jt C CM*, HIE 

• CCPPsASM-3 3. 

ipu.ec.;) i;, i3 

13 CCFF«CCPP+:SL*CFLCAT (I 1-1.) 

12 CALL C AG IM»Z»CCF*) 

1«.‘ FA°A = FRMCO (M ,Z ) 

JO C i v, 2 ) 

MR“TC/t!,12"j APS.PllOA 

12: FC^AT (ICxI* /•'PL. “ D T . >VTq“E s* ,r 11. 4 ,JSPU IS . TOT. £ NT AS a 

1F11.4) 

Xl=f'LOA T ( ? M 

F?sAL0GiC(*t)/UCGl3 13.3)+) .1 
PF-TM’ KM) 

1C" 3 

3 0 72 FC-lfAT SS (2 (lx, <<*£ 12. 3 ) >) 

CALL AFFTfZ, 

I C" '• 1 

CALL ArFTIZt I f> , 3 ,MM , 1 1 ) 

I'OU 

WAITstc, ,1:13) NC,Pa«»A 

iji’’ fop.va“ (i :x,<cclp;es ncs*,i4, ****cl:ss.uce totals entfse p 

1F11..A, /3cM/) 

IF< IT0.A4 tec. 1) 3.5,31 

30 B«:*)T 22 

22 FC^W (1CX ,*ENVELCF9® COFPLEXg S I a N A L. C,6NTR5g*> 

CALL TCACcNV It', 2, FACT, J> , *' > 

31 ITCTRAl *EC. 1) 63,7'* 


b5irs(e»i:) hCroacA 

FC«! W AT (i CO "PL CxE SICmi enr-ir „ ,. - _, 

1 *H.U S si? ssssc, «ous c ?un:lK?;^1 0 frj«,5'“*;;;- sy>) 

Cl ,3 .SE UE t »t CNAN t ILL Ch OE SIGNAL 

call c cu i f n ,z ,?/c,x » y , r j 

passage cans s/p a^pli pi 

cc/irms 

CC 123 L»i,»l 

cu*vv 

CALL Al'FKh, 2.w, A ,B» 

CC 124 L«1,N 
ULi'Z CL I* VO 

IPdrqAS.ES. 1) 32,33 
P“I NT 21 

SMMi SW!e ” 

IFttT»A8.5«.l, ,1,91 
PPINT 21, NC, PAPA 

FCRMAT (1 25! ,PSNVELCP°£ COPPt, £Xg SIGNAL ^OPTIc Kf=i tl 
TRACE SL* 3E.SSCN cql; "L'lS. TCT re * Pu , i ^ * 

CALL CCUSOm.Z.NC.X^T) * KT * c £ ^Pll.^tACQFje) 

CAUCLL TL SPuCTPE CE SC^TIE 

CC,*| T IMUf; 

IC" 1 

rri Trifr T i 7,i,iV >S * 

IP ( IT*t AS , J ) 34,35 jy; 

horrsi g, 2g) nc w 

FC’^ J Ar <1H1,1CX, XSPECT^E CE SCPT-c Nr-'- ru 

-r( ITRA3 *BC» 1 ) 13C,2JC 

T^ACs CU SPECT2E CE SCrTIE 

CCf! r lNLE 
^ « N T 111, NO 

r -’T«') .3 li ^ ” ' iP IC ‘ SOfii ’ i * = ^ , 14, fTPACE 8ENSCN») 

C'P» J .7 

C C !< 7 IN L E ’ ° 7 ' * E ( * Z * N ° ’ > ' V ^ » 0 p » I ■ 3 ? i , P T , 5 . F ) 

C Cl T IN li 

, ,EI ‘ :) G0 rc 3116 

mi . 3 

t;c m kai t j M <j 
^K" IM/I' to 
JK»KtUf» ( K> 

VVB(;4Ak-fi.vj h , /STSp * J# j 
K k*I a IP X ( r Y) 

corps.) , 
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u 


tz 

l 7 


C 

125 
125 , 


ir 


CC' 


i <i t 


CC I* 1 ,KKK 
CC 15 i,*!. 1 ! 

ZUJeZZt*} 

ccr.’“,*ar: 

cc ee iA t 'rv-?' 5 . 

.fj. i) 4<‘it’3 

c :p f »c c r * ♦ xt 5 1* * c ho a t < : w, » 

CALL CAO <H,Z rCCFP) 

ptpuck»(. r : m,z) 

/r/,«^A^;rN| z> 
pi: m spa»a 

CAL.CU. Kf-vaCFf»£ SZGhn ENTICE 
10*7 

CALl AfffTt2i2f.VtSiKM,IC) 

XC»-t 

CAl.l, A FF t ( & , X N V » ,"u i'H » I O 
hC'U'ALlSA - :')* 

CC 125 

ZtL}*Z (L ) * VV 

CALL A 0PZ(.S»2»I',A»0) 

CC 126 , L* 1 ,M 
Z (U sZ IU**VV 
IC'3 i 

CALL AFf T ( 2 r IftV 1 3 t I'M 1 1 O 
GT.SZ l/,<) 

6 C »2 CJM 
CCs-K/C-Z 


CCf 


• « u 


{ f.C ) 


CCU R 1 ( I> c« -n T »5 C* CO* CC £ # COO } 

C S V " 1 ts ) ; , * A L C '< 1 0 ( C s L I ( I ) I 

f J VJ|I , »> k « 

: B r*c:»c:s,-.\) cc T o izz 


* A » « * 

*«» I •* V # 


C C u. 

•’C . 
cc sf: 

irC:, cc re 6 5 
>iffS6uizi:sn 
YU A 7 .* AG • l izr } ) 

5CTC Ts rul TO 7 + /M *• v f> 

CCMTlfr.s 

pc** » i3 ,*41001: no tot > 

^<* 3 6‘ALt2(f,Sn 
Y*S 6CMAG (z (NCJ) 

*x *»* Y'< *YK 
,*5LC01?{SOK) 

ecu I;; n *»c>:-»po7 

cc iti.ilh 

FK'i- (1 :/ »* 7P1CS C'J AiPf-C-IT 
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